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Abstract 

The differential cross sections for low-energy muonic hydrogen atom scattering from hydrogenic 
molecules are directly expressed by the corresponding amplitudes for muonic atom scattering from 
hydrogen-isotope nuclei. The energy and angular dependence of these three-body amplitudes is 
thus taken naturally into account in scattering from molecules, without involving any pseudopo- 
tentials. Effects of the internal motion of nuclei inside the target molecules are included for every 
initial rotational-vibrational state. These effects are very significant as the considered three-body 
amplitudes often vary strongly within the energy interval ^0.1 eV. The differential cross sections, 
calculated using the presented method, have been successfully used for planning and interpret- 
ing many experiments in low-energy muon physics. Studies of nuclear capture in pfi and the 
measurement of the Lamb shift in atoms created in H2 gaseous targets are recent examples. 
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I. INTRODUCTION 



A calculation of the differential cross sections for low-energy scattering of muonic hy- 
drogen atoms from hydrogen-isotope (hydrogenic) molecules is the main subject of this 
paper. The cross sections are expressed in terms of the corresponding amplitudes for 
muonic atom scattering from hydrogen-isotope nuclei. Thus, a dependence of these three- 
body scattering amplitudes on the collision energy, scattering angle, and spin is directly 
included. For numerical calculations, the three-body amplitudes computed using the adia- 
batic method Q i, 3 0, H, Q H , are employed. 

Many experiments in low-energy muon physics are performed using molecular hydrogen- 
isotope targets (see, e.g., Refs. 3 planning and interpreting such 



experiments, the differential cross sections for the following processes are often required: 

elastic scattering: afi{F) + BC afi{F) + BC, (la) 
isotopic exchange: a/i + BC ^ bfi + AC, (lb) 

and spin-flip: a/i(F) + AB ^ a^{F') + AB. (Ic) 

A monic hydrogen-isotope atom in the IS* state is denoted here by or 6yu; F and F' 
stand for the initial and final total spin of the muonic atom. The molecules BC, AC, and 
AB denote the hydrogenic molecules H2, D2, T2, HD, HT, or DT. The processes (0) can 
take place with simultaneous rotational and vibrational transitions in a target molecule. 
Thus, the name "elastic" assigned here to the scattering (jlap refers solely to the state of the 
muonic atom. The cross sections for the processes (0) are henceforth called the "molecular" 
cross sections. 

For many years, only the cross sections for muonic hydrogen atom scattering from 
hydrogen-isotope nuclei ("nuclear" cross sections) were available. 

elastic scattering: o,fi{F) + b ^ aji{F) + h, (2a) 
isotopic exchange: afi + b ^ bfi + a, (2b) 

and spin-flip: a^{F) + a ^ a^{F') + a. (2c) 

The application of the nuclear cross sections to a description of experiments performed in 
molecular targets gives very unsatisfactory results. A characteristic kinetic energy of muonic 
atoms in typical gaseous targets is lower than a few eV ll^]. Therefore, it is necessary to 
take into account effects of molecular binding and electron screening. 

Since a muonic hydrogen atom is a small neutral system, the methods developed for 
the description of neutron scattering in matter can be adapted, to a certain extent, for the 
muonic atom case. Molecular effects in low-energy neutron scattering from nuclei bound 
in chemical compounds are estimated using the Fermi pseudopotential jl3, E, Such 
a pseudopotential is proportional to the constant scattering length. The Fermi method was 
used for decades for the calculation of low-energy neutron cross sections (see e.g., Ref. 
and references therein). In particular, a quantum-mechanical treatment of slow neutron 
scattering from molecular hydrogen and deuterium was presented by Young and Koppel [2l| . 

A method of calculating bin ding effects in the molecular processes (P), based on the Fermi 
approach, was derived in Refs. 12 2^- In particular, specific spin-dependent pseudopoten- 



tials were introduced for a description of a muonic atom interaction with a single nucleus. 
However, this method has a limited applicability since the nuclear processes (j21) involve sev- 
eral partial scattering waves [^, |^ even at low (~ 1 eV) energies, in contrast to low-energy 



2 



neutron scattering. Moreover, muonic atom scattering often changes strongly (e.g., pfi + p 
and t/i + t) with energy in the intervals comparable with the rotational thresholds of hy- 
drogenic molecules. A solution to this problem is to base a calculation of the molecular 
cross sections on the full nuclear scattering amplitudes, which include all the angular and 
energy dependence. The effective radius of interaction between a muonic atom and a nucleus 
is much smaller than the internuclear distance in a hydrogen molecule jlj. Therefore, the 
amplitude for scattering from two bound nuclei can be well approximated by a sum of the 
two corresponding amplitudes for scattering from isolated nuclei. In such an approach, it 
is necessary to take into account the internal motion of the nuclei inside a target molecule. 
This motion can be neglected for a molecule consisting of heavy nuclei. However, we are 
dealing with the lightest molecules and, therefore, the kinetic energy of nuclear motion due 
to zero-point vibration is on the order of 0.1 eV. 

In Sec. ini the amplitudes for the molecular processes ([Q) are expressed in terms of the 
amplitudes of the three-body reactions (0). The derived formulas depend on the momenta of 
internal motion of the nuclei in a target molecule. The differential cross sections for scattering 
from molecules are obtained in Sec. Hill using a harmonic model of molecular vibrations. 
Also, electron-screening corrections to the cross sections are given in this section. Some 
typical examples of the computed differential cross sections are shown in Sec. IIVI 



II. AMPLITUDES FOR SCATTERING FROM MOLECULES 

Let us consider a/i scattering from a molecule BC consisting of hydrogen-isotope nuclei 
b and c and two electrons. First, we assume that the nuclei b and c are different from the 
nucleus a, so that the scattering is spin-independent Also, electron screening effects 

are neglected in this section. The scattering lengths of the processes a/i + 6 (c) 0| are 
much smaller than the molecular diameter Rq ^ 300 (a^ denotes the Bohr radius of the 
a/i atom). The intexaction of a muonic hydrogen atom with nucleus b (or c) is important 
at distances -C -Rq Ql- Hence, it is assumed that a/x interacts with a single nucleus during 
the collision with the molecule. We also assume that the molecular bond is unperturbed at 
the moment of collision. Therefore, a/x collision with 6(c) is treated here as if this nucleus 
were free, except for its momentum distribution due to the molecular binding ji^]- This 
means that the amplitude for a/i scattering from a bound nucleus is the same as that for an 
identical free nucleus, provided the momentum of the relative motion is not changed. 

At large distances between a/i and BC, the initial ipo and final ipn coordinate wave 
functions of the system are as follows: 

V'o(r, r^, R-) = 0i(r^) ^o(R') exp(iko ■ r) , 
ipnir, r^, R) = 0/(r^) <?n(R) exp(?k„ ■ r) , 

where ko and k„ are the initial and final momenta of a/i; and <P„ are the wave functions 
of the initial and final rotational-vibrational states of the molecule BC. The corresponding 
wave functions of the IS* muonic atom are denoted by 0j and 0j, where the indices i and / 
refer to the processes (j21) with the nuclear scattering amplitudes /«/ jE0,SS0,0- In Fig- III 
the relative coordinates used for providing a description of the system are shown. The 
vector R connects nucleus b with nucleus c; r denotes the a/i position with respect to the 
center of mass (cm.) of BC; is the a/i internal vector. The vector stands for a/i 
position relative to nucleus b. 
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FIG. 1: (Color online) Relative coordinates used for the description of muonic atom a/i scattering 
from a molecule BC. The cross denotes a position of the center of mass of this system. 



It is convenient to express the amplitude for afi scattering from BC in terms of the 
momenta kg and k„ in the center-of-mass system (c.m.s.) for afi + BC. On the other hand, 
the calculated amplitudes for afi scattering from nuclei are functions of the initial and 
final p'f, momenta in the c.m.s. of a/i + b. Therefore, further evaluation of the molecular 
amplitudes involves the investigation of a transition between the "nuclear" and "molecular" 
momenta. First, we assume that both the nuclear and molecular scattering can be described 
in the Born approximation. The amplitude for afi scattering from b bound in BC is thus 
given by the following formula (in muonic atomic units e = h = fia^j, = 1): 

^TT J (4) 

X C(R) <P}{r,) V^'\n, r,) <p,{r,) <fo(R) , 

where is the reduced mass of the afi + BC system (the masses of the electrons are 
neglected) and f^a^l is the reduced mass of the a^i atom 

Ma^ = Ma + m^, M^oi = Mfc + M, . 

The potential of a/i interaction with a free nucleus b is denoted by V^^\ Using the relation 

r, = r + /5fcR , (3, = MJ{M^ + , (5) 
in Eq. ^ leads to the following factorization: 

^Sl = - ^ / exp [.(ko - k.) ■ n]v;:^\n) 

^ (6) 
X / d^R exp [i/3fe(ko - k„) ■ R] C(R)^o(R) , 

in which vjj^ denotes the "nuclear" matrix element 

V^fiv,) ^ I d\<P}{r,)V^'\n,r,)Ur,). (7) 

The first integral in Eq. © is the Born amplitude for a/i scattering from a free nucleus b 
times the factor A4/ fib, where fib stands for the reduced mass of the system a/i + b: 
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The second integral in Eq. ® is a form factor describing the binding of b in BC. 

In order to investigate a dependence of the molecular amplitude J-'q'!^ on the internal 
motion of b inside the molecule, a momentum representation of the wave function <P„ is 
introduced 

1 



R 



(27r)3/2 



d^Kn exp{iPbKn " R) gn\l^n), 

J d^R exp(-i/?feK„ ■ R) <P„(R). 



(8) 



(27r)3/2 

The vector k;„ is the momentum of the internal nuclear motion in the final rotational- 
vibrational state n. The analogous equations can be written down for the initial molecular 
state (pQ with the internal nuclear momentum Kq. Upon, substituting Eq. (jHl) into Eq. (jlj) 
one obtains: 

J d^r d^Rd^Kn d^KQ exp [i(ko — k„) • r] 



''On 



X exp [i(3b{Ko - k;„) ■ R] Vy{rb) 



Then, using new variables and 



l^'b , r, Mabc 

M-b + P,^Qb: 



^ Ma Mabc 

R = — Qb 



Mab 

Mabc = Ma + Mb + M,, Mab = Ma + Mb , 

in Eq. © and performing integration over the vector x = {PbMabc/Mc)Qb, one has 



(9) 



(10) 



-jrib) 



M 



d^Kn d^KQ / d'^Vb exp 



. f^b . , n l^b 



■ n 



X V}j\rb) exp 
X 6{ko — k 



+ l^n) gn^^l^n) go\l^o) 



'111 



The integral over r;, times —fib/^ir is the Born amplitude f^j^ for afi scattering on a free 
nucleus b, expressed by the momenta k„, k;„, ko, and kq. Thus, Eq. (|TT|l can be written 
down in the following form: 



[b) 



On 



f j3 ^3 Ab) f ^b . , n f^b A^b , . n 

J d Knd Kq tlf I — ko + /^/ft— HQ,—kn + (3b 



1 M 

„ .^^^^^ , -"Mfc-'-'A^"" ■ '^"M, '7 (12) 

X (5(ko - k„ - Ko + l^n) gn^*{l^n) fl-J^Vo) , 

where subscripts i and / label the kind of nuclear process (0). Now, we make the basic 
assumption that Eq. (fT^ is fulfilled by the exact nuclear amplitudes flj\ The integration 
over Kn is performed readily using the conservation of the total momentum, which gives: 



^0^? (ko, q) = 4 — / d'no flfiPb , Pb + q) g^^*i'^o + q) ^'^(^o) , (13) 

Pb f^b J 
l^b -, , r, l^b 



p,^-ko + /3,-.o. 



Pb = Pb + q 



(14) 
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The vector q denotes the momentum transfer 



q = k„ - ko = Pi - Pf, . (15) 

Thus, the molecular amplitude J-'^^J for scattering with a fixed momentum transfer q is 
determined by the free nuclear amplitude /-^^ with the same momentum transfer. However, 
the initial momentum Pf, in the afi + b c.m.s. is different from the initial momentum ko in the 
molecular c.m.s.. According to Eq (jl4|) . the vector p^ depends also on the internal motion 
of b. This gives the following afi kinetic energy in the afi + b c.m.s.: 

e, = -^e + P!-^ e,, + 2/?,^^ cos 9 , (16) 

where 6 is the angle between the vectors ko and kq. The muonic atom kinetic energies in 
the nuclear (e^) and molecular (e) c.m.s. are 

The internal kinetic energy of the molecule BC is 

ebc=:^, f^bc = M,-' + M-\ (18) 

At £ — > 0, the collision energy Sb in the a/i+b system is determined solely by etc- This 
energy never vanishes because of the zero-point vibration of the molecule. In particular, 
for the lightest H2 molecule, the second term of Eq. (fTB|) is on the order of 0.01 eV. This 
energy is inaccessible, but it affects the molecular scattering amplitude (fT^ . For a fixed e, 
the spectrum of Eb is quite wide. Its width is determined by the term y/esbc which depends 
on the hydrogenic-molecule vibrational quantum 0.3-0.5 eV). Therefore, at a given e, 
the molecular amplitude (jl3|) contains contributions from the nuclear amplitude fif taken 
at different energies. This effect should be taken into account when fif changes significantly 
within the spectrum ()16|) of Sb, which often occurs in muonic atom scattering. 




FIG. 2: (Color online) Relations between the initial and final momenta and between the scattering 
angles -db and t9 in the nuclear and molecular c.m.s. (for ko ^ ko). 



If ko ^ kq and kn ^ the dependence of on in Eq. (fn|) can be neglected. As 
a result, we get: 

P.^;^ko (19) 
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and Eq. ()13|) is then factorized as follows: 



•J n-, 



1 M 



On o-i 



flf (Pt , P6 + q) / d'Ko 9l!'> (Ko + q) g!'^ M. (20) 



This formula can also be used when /-^^ weakly depends on the variation of ph due to the 
characteristic spectrum of kq. After the substitution of the Fourier transforms (jH} into 
Eq. and integration over kq, we obtain 

-(b) _M Jb) f fJb 



X / d^R <?;(R) exp(i/3bq ■ R) . 



The molecular scattering amplitudes ()T3j) and (j2T|l naturally take into account the depen- 
dence of the nuclear amplitude fif on the a/i energy and scattering angle i). At low collision 
energies, the amplitudes fif are well approximated by the corresponding constant scattering 
lengths A°j. As a result, Eq. (j7H) is simplified: 

^o^^(q) ^ -— I ^n(i^) ^M^Pb^ ■ R) ^o(R), 

(22) 



M 



On iqj „ * ^ ^j/ • 



This equation can be formally obtained in the first Born approximation, using the pseu- 
dopotential 

* (r,) = 5L aIJ)" S(r + AR) (23) 



in Eq. Q. Such a potential for a constant scattering length was first introduced by Fermi [1 
and more rigorously derived by Breit and then by Lippmann and Schwinger 19l |. 

In the case of a general spherical potential with a finite range, it is possible to generalize 
the Fermi pseudopotential by the introduction of partial pseudopotentials corresponding 
to subsequent scattering waves ji^ 2^|. However, the calculations of the cross sections 
for muonic atom scattering on molecules, presented in this paper, are directly based on 
a knowledge of the amplitudes for nuclear scattering. Therefore, a formulation of a gen- 
eralized pseudopotential, in this case, is superfluous. Let us only note that the correct 
nuclear partial amplitudes can be formally obtained by the substitution of the following 
pseudopotentials: 

V}p-\n) = — (2J + 1) Xfl\pb) 5(r,) Pj(cos^9,) (24) 
l^b 

into Eq. (@]). The energy-dependent nuclear scattering length Xf^^ is defined as follows: 

f^f = -(2J+l)\f^'Fj{cosdb), (25) 

where /^^''"^ are the partial nuclear amplitudes for afi + b scattering 0, B IE S 0] and J is 
the angular momentum of the an + b system. The angle between the vectors and p'f, is 
denoted here by -^b- The function Pj is the Jth Legendre polynomial. 
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The problem of the angular and energy dependence of the nuclear scattering t/i + d in the 
t/i+D2 process was alternatively solved in Et'I by the introduction of the effective polarization 
potential (~ r~*^). However, the magnitude of such a potential for the tfi-d interaction was 
determined separately for every given collision energy. Thus, such an approach is more 
complicated than the direct use of the computed nuclear-scattering amplitudes and neglects 
a wide distribution of the deuteron kinetic energies in the D2 molecule. 

When the internal motion of the nuclei inside the target molecule cannot be neglected, 
the molecular amplitude is given by Eq. ()13p . However, in a general case, the numerical 
evaluation of the integrals over kq is difficult. The role of the internal motion is most 
important if the condition g <^ is fulfilled, which implies that the internal state of the 
molecule is not changed {gn^ = gQ^)- In this case, Eq. (fT^ is approximated as follows: 

M 



7 (26) 

with Vll'^ (kq) being a distribution of the momentum of nucleus b in the molecule: 

vl,'\Ko)^f3^'\gS'\^o)\'. (27) 

For elastic scattering at /eq — ^ 0, when one expects that the internal-motion effect is the 
strongest, the molecular scattering amplitude is 

^S' « -- ^, (28) 

where the bar denotes averaging over k,q. Equations and ()28|1 suggest a reasonable 
approximation of formula for finite ko. When the exact nuclear amplitude f^j^ is replaced 
by the averaged function 



f^^\ko, q) = - X](2J + 1) \f^-\ko) Pj(cos^,) , 

(29) 



'if 



Eq. (jl3p is factorized. This leads, finally, to an equation similar to Eq. (j2ip . with f-'r' replaced 



by the mean amplitude f-j\ This approximation gives the limit at ko — ^ 0. On the 
other hand, this approximation coincides with the asymptotic amplitude (pT|l at ko 3> kq- 
The dependence of {}b on Kq is neglected here since the higher partial waves (J > 0) in the 
nuclear scattering are important only at k^^ Kq. 



III. MOLECULAR DIFFERENTIAL CROSS SECTIONS 
A. Spin-independent scattering 

In the presented approach, the total amphtude J-'on for a/i scattering on a molecule BC 
is equal to the sum of the amplitudes for scattering on the bound nuclei b and c 

Toni^O, q) = 4ni^o, q) + 4ni^o, q) , (30) 



8 



where J^^'^ is given by Eq. (j^H), and the derivation of JFg^^ is analogous. Let us first consider 
the spin-independent case a ^ b, c. Assuming that vibrations of the molecule are harmonic 
and that there is no coupling between the vibrational and rotational degrees of freedom, the 
molecular wave function ^n(R-) takes the form: 

<P^{R) = ^YkmAR), ^ = (31) 

where quantum numbers K, Mk label the rotational state of BC. The radial wave func- 
tion corresponding to vibrational quantum number u is 



,{R) = K [a(i? - Ro)] exp [-^a^iR - Rq)^] , 

I 5 ^ (32) 



where H,^ denotes the z/th Hermite polynomial. The rotational and vibrational energy 
levels are given as 

Ek = B,,,K{K + 1), E,= {u+l)u;o. (33) 

At the temperatures usually applicable to experiments, hydrogenic molecules are initially in 
the ground vibrational state z/ = 0. 

Inserting the expansion of the free-wave function (in terms of the spherical Bessel func- 
tions j; and the spherical harmonics Yim) into Eq. (PT|) . one obtains for the bound nucleus b 



=47r — Yl ^^^^(^^1) ^irniq) 

l,m (34) 



df2R Y\,^,^{R) Yim{R) YkmAR) 
The real function V^i{Phq) is a result of the integration over R 

POO 

VMq) = I dRu^iR) iiiPtqR) uoiR) , 
V,i{(3bq) 



1 if z/ = and / = ^^^^ 



jo otherwise. 

The initial state of the molecule is denoted here by the set of rotational and vibrational quan- 
tum numbers 0={K, M^, z/=0). The final state is labeled by n={K', M^, u). The indices i 
and / in the nuclear amplitude /^^^ refer to the initial and final states (spin or isotopic) of 
the scattered muonic atom. Integration of the three spherical harmonics in Eq. ()34|) over 
the solid angle leads to the following result: 



= [47r{2K' + 1){2K + 1)] (-1)*^^ t^-^ 
M 



K' I K\ ( K' IK 
j \-M'k m Mk 
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expressed by the Wigner 3j symbols. 

The molecular differential cross section, averaged over the projection Mk of the initial 
angular momentum and summed over the projection M'j^ of the final angular momentum, 
is equal to 

daon 



2K+1 ^ I ' 



mol I 2 
On 



(37) 



Mk,M'k 



The solid angle fl{d, ip) is connected with the direction of the vector k„ with respect to the 
initial a/x momentum ko (see Fig. 12)). Substitution of Eq. (jHUj) and the analogous formula 
for the nucleus c into Eqs. (jHUj) and (jHTf) gives the following cross section: 



On 



dQ ko 



K'lK 



where the angular-momentum factor Wk'ik is defined as 



(38) 



Wk'ik = {2K' +I)i2l + l] 



The reduced mass of the a/x + c system is 



K' I K 




Km 



and 



When BC is symmetric (6 = c, fib 
Eq. takes the simpler form: 



-f3b = Mb/{Mb + M,). 

fic = fJ', Pb = Pc = P 



2' 



and /; 



if 



(39) 



Ac) 

J if 



dai 



On 



dn 



<^)V.iti:[M)-+i] 



(40) 



xWK'iKVliiPq) 



The molecular cross sections (|38p and (j40|) directly include the dependence of the "bare" 
nuclear amplitudes on the collision energy Sb (^c) and on the scattering angle 'db (t?c)- 
They are derived for high collision energies kq. However, they can also be used at lower 
energies as a reasonable approximation if the nuclear amplitudes ff^^ ifif) replaced by 

the amplitudes fl^^ ifif ) averaged over V^'^ (V^o ^). In the case of the wave function (jHlj) . 
the momentum distribution (j27|) for the ground vibrational state i>=0 of the molecule BC 
has the form 



Pb -^0 -2 



j;^(/3feKoi?o)exp 



Pb'^O 



(41) 



X Y 



KMk 
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in which kq = kq/kq. After averaging Vq^"^ over orientations of the molecule, one obtains 



Ji^lPb^^o-KoJ exp 



3/2 



(42) 



A distribution Vqk of the internal kinetic energy of the target molecule can be derived 
similarly. For z/ = 0, one has 



VoK{£bc) dSbc = — z — j/cl-Roav^) exp(-t^bc) 



vr 



(43) 



be ■ 



where Uhc = Se^c/'^o- This distribution is widest for the lightest H2 molecule. According 
to Eq. (|TT)jl . this leads to a broad distribution of the collision energy Sb in the nuclear c.m.s. 
for a fixed collision energy e in the molecular c.m.s.. In Figs. El and EJ the calculated 
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FIG. 3: (Color online) Distribution of collision energy in the pfi + p c.m.s., for a fixed collision 
energy e in the pfj,+}i2{K = 0) c.m.s.. 
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FIG. 4: (Color online) The same as in Fig. El for e = I and 10 eV. 
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Eb spectrum for the ground-state H2 molecule, at several values of e is presented. At the 
lowest energies, the shape of this spectrum is mainly determined by the motion of a bound 
proton. The mean value of Eb equals 0.036 eV for e = 0.001 eV. For e = 0.1 eV, a single broad 
peak with a mean value of 0.11 eV is observed in the Eb spectrum. Only at e > 1 eV, does the 
average Eb approach the asymptotic value of {jibl M.)e given by the first term of Eq. (fTB|). 
However, the width of the Eb distribution, determined by the last term of this equation, 
increases with rising e. The ratio of this width to the mean value of Eb decreases as 



B. Spin-dependent scattering 

When at least one of the nuclei bound in BC (e.g., the nucleus h) is identical with the 
nucleus a, it is necessary to consider spin-dependent reactions (^aj) and (Pc|) . Let us introduce 
the following notation 



J" = + Sa + Sb + Sc , 

(44) 



F = + , I = S5 + Sc , 



Sb = F + Sb , S, = F + Se , 

= F + I = Sfc + s, = Sb + Sc , 

where s^, s^, s^, and Sc, are the spins of the muon and of the nuclei a, 6, and c, respectively. 
Thus, JT" is the total spin of a/i + BC system and I is the total nuclear spin of the molecule. 
It is assumed that is conserved in the scattering since the spin-orbit interaction is very 
weak. Also, it is assumed that the spin Sf, (Sc) is conserved in local collisions of a/i with the 
nucleus h (c) bound in BC . 

If the isotope c is different from a and 6, the molecule BC is asymmetric and its parity 
is not definite. Thus, the directions of the nuclear spins Sf, and Sc are independent of each 
other, and a unique spin S;, of the subsystem aii + b is assigned to the initial and final ^/ 
spin states of the system a/i + BC (with fixed values of F and F'). These states can be 
written down as follows: 

e.(Sfe,s,;F)=eaM_fe(Sb;F)ec(Sc), , 
^f{Sb,s,;F')=^a>.-b{Sb;F')U^,), ^ ^ 

where ^a/^-b and are the eigenfunctions of the conserved spins Sb and Sc, respectively. In 
this case, the total wave functions take the form 

■ipo = 0i(r^) ^o(R) exp(iko ■ r) ^^(8;,, s^; F) , 
V^n = 0/(r^) ^n(R) exp(2k„ ■ r) ^/(Sfe, Sc; F') , 

at large distances r. In order to obtain the correct molecular-scattering amplitudes, sym- 
metrization of the functions ipo, "^n over the two identical particles a and b should be per- 
formed. As a result, we obtain the molecular amplitude JFq^-* which is expressed by Eq. flHH|) 
with f^^^ replaced by the spin- dependent nuclear amplitude fpp, for the process (^aj) or 
The calculated amphtudes fpp, are already symmetrized over the identical nuclei 0, 0, 0, 
Similarly, the total molecular cross section c is given by Eq. (j38p 

with f-^^ replaced by fpp,- In the case of spin-flip reaction, we substitute f-^^ = in Eq. 
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because this process is very weak when the isotope b is different from a [28?]. The cross sec- 
tion da^^/dQ can be averaged over the projections of spin Sb, which gives the mean cross 
section daon/dQ. 

When all the hydrogen isotopes are identical, the molecule BC is symmetric and its initial 
and final parities P/ and P[ are definite. Therefore, the values I and /' of the molecular 
spin are definite. The total spin functions of a/i + BC are now eigenstates of with fixed 
values of F and / (or F' and /'). Thus, the total spin state is determined by four quantum 
numbers: the absolute value J of the total spin J^, its projection J7^, F, and /. The initial 
and final total wave functions of the system with a = b = c are 

V'o = 0i(r^) M^) exp(iko • r) ^J; F, I) , 
^„ = 0/(r^)^„(R) exp(^k„-r)0(J^;i^',/'), 
with the following condition to be satisfied 

F + I = F' + I' = :7'. (48) 

The functions and ^/ contain contributions from different states of the operator Sf,. They 
can be expanded as follows: 

CiiJ; F,I) = J2 CbiSt, se; F, I) Ca,-b{Sb) U^c) , 

X (49) 
^fiJ; F', I') = J2 C'biSb, s,; F', I') ea^_,(S,) U^,) , 

Sf, 

with Sc subject to the condition: S^+Sc = J'. The factors Cb{Sb, s^, F, I) and C^(Sb, s^, F', I') 
are obtained by the multiple use of the Clebsh-Gordan coefficients. The expansions of the 
total spin functions in terms of ^afi-bi^b) are necessary since the presented method is based 
on knowledge of the three-body scattering amplitudes fpp, evaluated for fixed values of Sb- 
After performing a symmetrization of the total potential V^'-''^ + V^'^^ and of the wave 
functions ()47|) over the three identical nuclei and proceeding as in Sec. |n] for the spinless 
case, one obtains 

^On = ^ Cb{Sb, Sc, F, I) Cl{Sb, Sc, F', I') fpp, 

^ (50) 



X j d^R$l(K)exp{iPq- R) ^0 (R) 



In the derivation of Eq. ()50|) . it has been assumed that the three nuclei are never close 
together, i.e., the nucleus c is only a distant spectator when a/i collides with the nuclei b. 
As a result, the molecular amplitude (j5(jp is expressed in terms of the three-body ampli- 
tudes fpp,- By employing Eq. (j^Uj) . the total molecular amplitude takes the form 

J'on =—dFF'{K,K') [ d'R<P:{K)[exp{tf3q-K) 

J (51) 



where ^ff'{I,I') is given as 

W(/, /') = Yl ^"^^"^ ^' ^) se; F', I') f%, . (52) 
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This result is independent of a choice of nucleus b because of the symmetry: b c, Sb Sc, 
and Sc ^ Sf,. Using the expansion of the plane wave in terms of the spherical harmonics 
in Eq. (j^Tj) and taking into account that [l + P/P/'(— 1)'] = 2 for every allowed rotational 
transition, we obtain 

in the case a = b = c. Since, in experiments, both the muonic atoms and the target molecules 
are not polarized, the cross section is averaged over J' and / and summed over /' (for fixed 
K and K'), which is denoted by the horizontal bar over the squared amplitude 



C. Electron screening corrections to molecular cross sections 



The differential molecular cross sections derived in the previous section include only the 
muonic-atom interaction with nuclei. It is necessary, however, to include electron screening 
effects in a/i scattering from molecules. At e < 1 eV, the relative velocity of afi and BC 
is smaller by several orders of magnitude than the muon velocity in afi and is also smaller 
than the electron velocity in the molecule. Therefore, it is possible to introduce an effec- 
tive electron-screening potential, which is obtained by averaging the Coulomb interaction 
between a/i and the electrons over the muon and the electron coordinates. The range of 
a dominant fraction of the afi-b potential is smaller than about 20 fi\. On the other 
hand, the a/i interaction with the electrons is important at distances on the order of the 
Bohr radius ~ 207 of the electronic hydrogen atom. Thus, a/i collision with an ordinary 
molecule can be described as scattering on the two potentials with very different ranges. 

The effective screening potential Vei for a/i scattering from hydrogenic molecules has the 
following form 2^, 2^: 



el 



exp ( r + /SfeR 



-f 25'^ exp 



-^{\y + (3M + |r-/3cR| 
a« V ' ' ' 



-Fexp ( - — |r-/3cR 



(54) 



in which 
and 



C = 2 + 8.4 



Sr,= {l+Wr, + \ vol) exp(-W,,), 



= {Ma-m^)/{Ma + m^), 
rjRo 



7] = 1.2 



The electronic correction to the process (fTa|) is calculated using the first Born approximation. 
The total molecular amplitude J^on°^ equal to sum of the nuclear amplitudes and the 

screening amplitude jFg^ 

^or(ko, q) = q) + ^o^^ (ko, q) + ^ot(q) , 



(55) 
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where 

•^ot(q) = / d'R exp(-zq ■ r) ^^^^ 

x<?;(R)v;i(r,R)^o(R). 

The calculated amplitude JFq^ falls rapidly when qa^ > 1, which occurs even for the lowest 
rotational excitations of a hydrogenic molecule. Therefore, it is sufficient to take this am- 
plitude into account only for the strictly elastic scattering. The isotopic exchange (j2E|l and 
the strong spin-flip (Pcj) reactions are due to the exchange of the muon between two nuclei 
taking part in direct collision. Therefore, a/x scattering from electrons cannot cause these 
reactions. The first non- vanishing screening corrections to the spin-fiip or isotopic-exchange 
cross sections appear only in the distorted wave Born approximation. 

A further evaluation of the screening corrections to the scattering amplitudes should 
be performed numerically. With regard to the elastic processes, these corrections are very 
significant. For example, at £ — > 0, the screening amplitude JF^^ for p/i + H2 elastic scattering 
is comparable to the corresponding pfi + p scattering amplitude. In this limit, the relative 
screening corrections to the spin-fiip or isotopic-exchange amplitudes are on the order of 10%. 
At e > 1 eV, screening effects practically vanish for all processes. 



IV. EXAMPLES OF MOLECULAR CROSS SECTIONS 

In this section some typical examples of the molecular cross sections are shown. The 
nuclear scattering amplitudes given in Refs. 0, 0, IE 0, 1^ are used as the input for 
computation of the molecular differential cross sections. These amplitudes are first averaged 
over the internal motion of nuclei inside the target molecules, according to Eq. ()29|). 

^ 1 
B 

o 

2 0.9 
o 

i 0.8 
0.7 



0.6 
0.5 

0.1 0.2 0.3 0.4 0.5 0.6 
energy e (eV) 

1/2 1/2 

FIG. 5: (Color online) Amplitude /^{ and average amplitude fi{ for = 0) + p scattering 

versus collision energy e in the molecular c.m.s.. 




1^ , ■ . , 1 /2 

In Figs. IHl and El the averaged nuclear amplitudes f^l for the elastic scattering = 

1 /2 

0)+Il2 and tfi{F = 0)+T2 are shown. The input nuclear amplitudes fi{ are plotted versus 
collision energy e in the molecular c.m.s., using the high-energy asymptotic relation e = 
{M./ lib) £b- The cusp in the elastic cross section pfj, + p, located at the spin-fiip threshold, 
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' I 

0.1 0.2 0.3 0.4 0.5 0.6 

energy e (eV) 



FIG. 6: (Color online) The same as in Fig. El for tfi{F = 0) + 1 scattering. 



is smeared out after the averaging over proton motion in H2. Although the vibrational 
quantum for T2 is smaller than that for H2, smoothing of the amplitude for t/i + t elastic 
scattering is also important, owing to strong changes of its value within the energy interval 
of 0.1 eV. This is particularly visible in the cusp region and at £ — > 0. This smearing strongly 
affects the molecular cross sections, which are expressed by the squared amplitudes. The 
elastic cross sections for dfi + d scattering are quite flat at the lowest energies 
a result, differences between the amplitudes /^^ and ffj are much smaller than in the protium 
or tritium case. The role of smearing effects were investigated during final a analysis of the 
PSI diffusion data A spectacular improvement of the fits to the data, especially for p/i 
diffusion in H2, was achieved when the averaged nuclear amplitudes ffj were used for the 
calculations of the molecular cross sections. This mainly concerns the elastic cross sections, 
as smoothing effects are generally much smaller in the case of spin-flip or isotopic-exchange 
amplitudes, which weakly depend on the energy below 1 eV 0, MSB 111- 

Electron-screening and molecular-binding effects are clearly seen in the molecular differ- 
ential cross sections. The range of the screening potential (j54j) is on the order of Cg, so 
that the condition k^ae ~ 1 is fulfilled already at £ ~ 0.001 eV. Many partial waves begin 
to contribute significantly to the screening amplitude above 0.001 eV. As a result, the 
molecular cross sections are anisotropic even at very low energies. Moreover, scattering 
from a molecule is connected with different rotational transitions, which additionally leads 
to a complicated angular distribution of the scattered atoms. This is in contrast to a/i scat- 
tering from a bare hydrogen nucleus, where few partial waves contribute significantly to the 
nuclear cross sections below 100 eV. In all of the three-body cases, the s-wave cross section 
describes scattering below 0.1-1 eV well, which is therefore isotropic in the nuclear c.m.s.. 
The cross section daQn/df2 for p/x scattering from a ground-state H2 molecule is presented 
in Fig. [7J for e = 0.25 eV and e = 0.5 eV. Scattering at the angles {} < 30° is dominated by 
pfi scattering from the electron cloud. The peaks at greater angles are due to the rotational 
transitions K = ^ K' = 2 and K = ^ K' = 4. The angular positions of the scattering 
peaks change with the variation of pfi energy. 

The differential cross sections for pn scattering from a free proton and from a H2 molecule 
are plotted in Figs. |H1 and El as functions of the collision energy e and the scattering angle 
in the molecular c.m.s.. Figure |H1 illustrates the cross section for the elastic pfi + p scattering, 
multiplied by 2. The mass of the target particle is, however, set to the H2 mass, for the 
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FIG. 7: (Color online) Differential cross section for pfi{F = 0) scattering from a H2(i^ = 0) 
molecule versus scattering angle -d. 




FIG. 8: (Color online) Doubled differential cross section for pfi{F = 0) + p scattering versus e 
and ■d. 



sake of comparison with this molecular target. Only the s wave contributes to the pjj, + p 
cross section in the considered energy interval. However, the scattering is not isotropic in 
the molecular c.m.s.. A cusp is apparent at the energy of the spin-flip threshold. The 
surface in Fig. El describes the corresponding cross section daon/df2 (cut above a value of 
3.5 X 10"^'' cm^/sr) for the ground-state H2 molecule. The electron-screening contribution 
to the differential cross section is clearly seen at all angles for e < 0.05 eV and only as 
a forward peak at higher energies. The lowest rotational transitions can be distinguished 
in this plot. Also a smearing of the molecular cross section due to the proton internal 
motion in H2 is visible, especially in the spin-flip threshold region. The sections of the 
surfaces shown in Figs. |H1 and El with the plane e = 1 eV are already quite similar, apart 
from the forward scattering. Thus, the molecular scattering at larger energies and angles, 
which contains contributions from many rotational- vibrational transitions, approaches the 
"doubled" nuclear scattering. For a better comparison, the angular dependence of the 
doubled nuclear and the molecular cross sections are shown in Figs. ^1 and ^2 at a fixed 



17 



FIG. 9: (Color online) Differential cross section for pn{F = 0)+Ii2{K = 0) scattering versus e 
and -i?. 
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FIG. 10: (Color online) Angular dependence of the differential cross sections for pfi{F = 0) scat- 
tering on a proton (doubled) and on a H2(i<' = 0), at e = 0.1 eV. 



collision energy. For e = 0.1 eV, the two cross sections are very different at all angles. For 
e = 1 eV, these cross sections are already quite similar at larger angles > 20°. One sees 
that there are only small rotational oscillations of the molecular curve around the doubled 
nuclear curve. 

The differential cross sections for the downwards spin-flip reaction in pfi{F = 1) scattering 
from a proton and from a H2 molecule are shown in Figs. ^1 and E| Electron-screening 
and molecular-binding effects in the spin-flip reactions are not so important as in the case of 
spin-conserving scattering (cf. Figs. Inland IT^. owing to higher momentum transfers. Larger 
differences between the nuclear and the molecular spin-flip cross sections appear mainly at 
small collision energies and small angles. 

Simulations of experiments performed in gaseous hydrogenic targets require knowledge of 
the differential cross sections for the molecular processes (^Q). These molecular cross sections 
have been computed and stored as computer files. They have been applied for planning and 
interpreting many experiments in H-D gaseous targets. For example, optimal conditions for 
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FIG. 11: (Color online) The same as in Fig. IIUI at e = 1 eV. 




FIG. 12: (Color online) Doubled differential cross section for the downwards spin- flip reaction in 
PIJ,{F = 1) scattering from a proton versus £ and i}. 



studies of n~ nuclear capture in p/x 11, 12| and for the measurement of the Lamb shift in pfj, 
atoms jisl . l3 . isj l created in H2 targets have been determined using the calculated molecular 
cross sections. These experiments are now underway at the Paul Scherrer Institute. 



V. CONCLUSIONS 

A method of calculating the differential cross sections for low-energy muonic atom scat- 
tering from hydrogenic molecules has been developed. This method directly uses the cor- 
responding amplitudes for muonic atom scattering from hydrogen-isotope nuclei, calculated 
within the framework of the adiabatic method for the three-body problem with the Coulomb 
interaction. Thus, the presented method naturally includes the angular and energy depen- 
dence of the three-body amplitudes in the scattering from hydrogenic molecules. Since, in 
many cases, the considered three-body scattering amplitudes depend strongly on the colli- 
sion energy within the interval < 0.1-eV, a broad distribution of the nucleus kinetic energy in 
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FIG. 13: (Color online) Differential cross section for the downwards spin-flip reaction in pfi{F = 1) 
scattering from a H2(i^ = 0) versus e and 



a hydrogenic molecule is taken into account. The molecular vibrations are described in the 
harmonic approximation. Therefore, the evaluated cross sections are valid below a few eV. 

For a fixed energy e of a muonic atom collision with a hydrogenic molecule, the calcu- 
lated collision energy Sb in the system consisting of the atom and a single hydrogen-isotope 
nucleus has a wide spectrum. At e —>■ 0, this spectrum reveals a shape of the kinetic energy 
distribution of the nucleus in a given rotational- vibrational state of the molecule. This effect 
is very significant, even for the molecular ground state, as the energy of zero-point vibration 
in hydrogenic molecules is quite large. The width of the ei, spectrum for the lightest H2 
molecule is on the order of 0.1 eV. At higher e, this width is even larger. As a result, the 
three-body amplitudes are strongly smoothed when calculating the molecular cross sections. 
This effect and the energy and angular dependence of the three-body amplitudes are included 
in the calculated set of the differential cross sections for low-energy scattering of IS* muonic 
hydrogen atoms from hydrogenic molecules. These are the only to date theoretical cross 
sections which give good agreement with many experiments involving and dfi scattering 
in gaseous H-D targets. The presented method can also be applied for scattering of other 
ground-state exotic atoms or neutrons from hydrogenic targets. 
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